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Abstract: 

Using the framework of QCD sum rules we predict the twist-3 contribu- 
tion to the second moment of the polarized nucleon structure function g2{x). 
As the relevant local operator depends explicitely on the gluon field, we em- 
ploy a recently studied interpolating nucleon current which contains three 
quark field and one gluon field operator. Despite the fact that our calcula- 
tion is based on the analysis of a completely different correlation function, 
our estimates are consitent with those of Balitsky, Braun and Kolesnichenko 
who used a three-quark current. 
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Considerable attention has been paid to measurements of the polarized 
structure function g\ in deep inelastic scattering. The unexpected EMC [[! 
results for the first moment of gf, the structure function of the proton, in- 
voked tremendous effort to gain more knowledge about polarized scattering. 
As a result of these intense discussions it became clear that polarization phe- 
nomena might offer much better opportunities to test QCD than unpolarized 
experiments. Meanwhile SMC || and SLAC || experiments have provided 
date on g™ and much improved ones on g\ . Mainly due to this experiments 
the general interest shifted to related questions most notably to the prob- 
lem of Q 2 dependence. It turned out that an estimate of the magnitude of 
higher twist contributions is urgently needed in order to compare the exper- 
iments of SLAC, SMC and EMC which cover different Q 2 regions ||. This 
is especially important for the SLAC data which correspond to Q 2 around 
2 GeV 2 . Besides being of great importance for the interpretation of the ex- 
periments it became also clear that the leading higher matrix elements test 
fundamental properties of the nucleon. Most notably on which we will 
focus in this contribution is a measure for the mean colour-magnetic field 
in the direction of a transverse nucleon spin. (Due to colour symmetry the 
expectation value (B a ) = 0, but the expectation value of the colour singlet 
(qBiq) = (q'joGoiq) ~ can be nonzero.) 

The first moment of gi(x, Q 2 ) at fixed Q 2 is given by || 

o ' dx 9l (x, Q 2 ) = i a <« + g (a® + 4d^ + 4/ ( 2 >) + . (1) 

The above formula does not include perturbative corrections (QCD). Higher 
order corrections were recently calculated in the leading twist approximation 
up to 0(a|) for non-singlet quantities, the Bjorken sum rule and - up to 
0(a 2 s ) - for the Ellis- Jaffe sum rule |J. The Ellis- Jaffe sum rule has a flavour 
singlet contribution. 

The reduced matrix elements for the twist-3 and twist-4 components of 
Eq. (fj) are defined by: 

1 (pS\q(0) h a g&' J + j P gG aa ] q(0)\pS) - traces = 



6 

= 2d® 



- (2p a p l3 S a + 2pPp a S a - p l3 p a S a - p a p a S p - p a p a S 13 - p a p /3 S 
6 v 

—traces] , (2) 



2 



and 

ip3\q(0)gG a ^q{0)\p3) = 2m 2 N f^S a (3) 

respectively, where \pS) represents the nucleon state of momentum p and 
spin S, S 2 = —m 2 N . 

The matrix element does not only lead to a contribution to the Ellis- 
Jaffe and Bjorken sum rules but is also in principle measureable as the twist-3 
contribution to the second moment of g 2 0: 

/ dxx 2 g 2 (x,Q 2 ) = -- f dxx 2 g 1 (x,Q 2 ) + -d {2) . (4) 
Jo 3 jo 3 

Although both twist-3 and twist-4 matrix elements are interesting, in 
the present paper we concentrate on the calculation of the twist-3 matrix 
element d^ using the QCD sum rules. Some time ago d^ and were 
estimated with this technique by Balitsky, Braun and Kolesnichenko [] (BBK) 
||. They got rather suprising results, namely a strong asymmetry between 

g?bbk (P r °t° n ) an d rfBBK( neu ^ rori ) an d a *^bbk (proton) which has opposite 
sign to that suggested by bag model calculations. ||. It has been argued 
[ |T0|j that these qualitative properties should be observable in other spin- 
dependent effects, like single spin asymmetries, such that their unambigious 
determination is of great importance. Any calculation of a non-perturbative 
quantity like d^ has to be based on a trustworthy well defmded method 
which is insensitive to technical details and physicaly irrelevant assumptions. 
In our opinion the only such available techniques are QCD sum rules and 
lattice gauge calculations. The latter are tried but are still in prelimenary 



state. [II]. Thus the reliability of the BBK result can only be checked 
by consistency tests, namley by comparing the results of calculations with 
substantially different interpolating currents for which the systematic errors 
are uncorrelated. The spread of obtained results gives a good estimate for the 
size of the systematic errors involved. There is actually reason for concern. 
The chiral structure of the three-quark current used by BBK is such, that 
graphs of a certain topology "accidentally" vanish because a trace of an odd 
number of gamma-matrices has to be taken. Unfortunately their sum rule is 
primarily determined by just a single graph and for this graph the up-quark 
contribution vanishes due to the special structure of the used three-quark 
current. Therefore the question arises if the obtained result reflects only 

1 In their notation d {2 ^ =<< V » /4 and / (2) = - << U » /m 2 N . 
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the the symmetry of the choosen current or if there is some deeper physical 
reason behind this. In this situation independent tests with different (more 
complicated) interpolating nucleon fields are obviously necessary. 

Unluckily the QCD sum rules technique used to estimate the matrix ele- 
ments present in Eqs. @ and (0) leads to very complex computations. As 
the consequence the fundamental results of BBK have not yet been checked. 
Given their importance for the phenomenology we have decided to meet the 
challenge and to perform an independent calculation with a different inter- 
polating current for the nucleon. 

The starting point for calculating the reduced matrix elements of an op- 
erator Or is the three-point correlation function 

n r (p) = i 2 j d 4 xe^ x J d 4 y{0\T{r]{x)O r {y)r}{0)} |0> (5) 

where r](x), (0\rj(0)\pS) = Xu(p,S) is an interpolating current for the nu- 
cleon. There exists a certain freedom of choice for this current. While all 
reasonable choices should give similar results the optimal choice must be tai- 
lored to the problem considered. For a long time the standard choice for 
sum rule calculations of nucleon properties has been the three-quark current 



introduced by Ioffe fT2 



rji(x) 



u a (x)C^u b (x) l5 Yd c (x)e abc . (6) 



This current is used by BBK in their calculation. 

As it can be seen from @ and (||), the operators defining and 
explicitely contain the gluonic field operator. This gluonic field has to be 
matched by another gluonic field operator. Because the three-quark current 
rji(x) does not contain the gluonic operator explicitely the desired contribu- 
tion must be generated through the perturbative emission of an additional 
gluon. The amplitude for this process is proportional to the strong coupling 
constant g. 

Alternatively, one may consider an interpolating field which in addition 
to three quark fields contains the gluon field explicitely. In that case gluon 
emission from the three-quark configuration has a non-perturbative charac- 
ter. 

A possible construction of an interpolating current with such properties 
was discussed in IO]. It was noticed that all the nice features of the current 
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in (|6|) are preserved by making it nonlocal, e.g. by shifting the <i-quark to 
the point x + e (and adding the required gauge factor to preserve gauge 
invariance). Expanding two times in e and averaging over the directions in 
Euclidean space e M e t 



\S^ u e 2 one arrives at 



Vg( x ) 



u a (a;)C7X^))757/^ Q /3 gG a(3 (x)d{x 



(7) 



Finally, projecting out the isospin | component leads to the following proton 
interpolating current: 



Vg(x) = 1(v^(x)-Vg( x ) 



where 
and 



_abc 



(u a (x)C ltM u\x)) 757^ [G a ?(x)d(x) 
(u a (x)C^d\x)) lf rf°<# [G^(x)u(x\ 



(8) 



(9) 



(10) 



We would like to stress here that the technique of QCD sum rules does 
not require the use of "the best" current from all the possible ones, it is only 
necessary that the current is not too bad in order that the contribution of 
interest is not suppressed by some special reason. 

In ref |H| the current @ was analysed in details and the corresponding 
overlap Xq defined through 



(0\ V (0)\pS) = m 2 N \ G u(p,S) 



(11) 



was calculated. Subsequently, it was used to predict the nucleon momentum 
fraction carried by gluons to be ~ 0.4 at /i 2 ~ 1 GeV 2 , and the gluon form 
factor of the proton at Q 2 ~ 2 - 3 GeV 2 . Both predictions agree very well with 
present day data and its generally accepted interpretation. This convinced 
us that the current t]g{,x) is an adequate choice to calculate properties of the 
nucleon. 

From the purely technical point of view the current rjc{x) differs in one 
important respect from the Ioffe current t]i(x). As mentioned before, the 
chiral structure of r)i(x) is such that graphs of certain topologies vanish be- 
cause a trace of an odd numbers of gamma matrices has to be taken. The 
same happens obviously also for r] old (x) but from correlators involving ?7 ox (x) 
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we obtain non-vanishing contributions form the graphs of the same topology. 
The calculation is therefore quite different from that described in ||. 

It is not yet possible to calculate the three-point correlator defined by 
Eq. (H) directly However, the usual philosophy of the QCD sum rules is to 
relate it to some other correlators using Wilson's operator product expansion 
(OPE). As it has been noticed in |14| for the case of a three-point correlation 



function the expansion is more complicated than in the case of two point 
functions and has the following structure: 

n r (p) = < o\ > +E^ L WnF L n (o) . (12) 

n n 

As usual the main objective is to separate different scales: vacuum expecta- 
tion values of local operators (0 L ) and correlators U BL describe long distance 
effects while the coefficients c L and c BL receive contribution only from highly 
virtual quark and gluon fields which propagate for small distances. The first 
part of the expansion - so called local power corrections (LPC) - is well 
known from two point QCD sum rules. The second part - bilocal power cor- 



rections (BPC) - has been introduced in [14[ to describe the large y region 
of the integrand in Eq. (El). This region cannot be properly described by 
local power corrections alone. Yet, the distinction between large and small 
distances is somewhat arbitrary and the corresponding scale is implicitly in- 
cluded in the definitions of local and bilocal power corrections. Changing the 
scale we can shift part of the BPCs to LPCs and vice versa. The optimal 
scale is always dictated by the nature of the physical phenomena one wants 
to investigate. 

The (OPE) of the correlation function ITr is given by the sum of both 
LPCs and BPCs, as in (|P2|). We note that in general only this sum has a 
physical meaning and is independent of the regularisation scheme. 

The bilocal power corrections are determined by the long-distance, non- 
perturbative contribution to correlation functions at zero momentum ||15|| : 

nj£(g) = i J d*ye^(0\T{O r (y)O n (0)} |0> (13) 

with q = 0. The Operators O n (0) are defined by the following series 

T(r,j(x)fj G (0)) = £ C BL (x 2 )O n (0), (14) 
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where the above expansion goes over a series of local, gauge-invariant op- 
erators of dimension n O n (0). When this expansion is inserted back in (|J) 
it gives the second contribution to Eq.flT^p. We have performed a detailed 
analysis of the BPC's to the sum rule for the twist-3 matrix element and 
have found that numerically important correlators do not contribute due to 
symmetry. In this respect the situation is entirely the same as in the case of 
the calculation discussed in ||. 

The QCD sum rules for the twist-3 part of g\ can be represented as 



= 2 d^XjX G S [tT p {lx] p u} m% + . 
(m 2 N ~p 2 ) 2 

[. . .] means antisymmetrization of the indices, while {. . .} denotes symmetriza- 
tion resp (see Eq. (||)). Here, we have singled out the double-pole due to 
the nucleon contribution in the intermediate state; other terms denoted by 
ellipses at the phenomenological right hand side are less singular in m 2 N — p 2 
and do not carry any useful information - later on they will be suppressed 
by suitable treatment. The couplings of the Ioffe current respectively our 
quark - gluon current to the the nucleon are denoted by Aj and Ac while 
denotes the nucleon mass. 

The left hand side has been expanded and all LPC contributions in- 
cluding operators up to dimension 10 have been calculated. It has required 
calculation of more than 100 Feynman diagrams. Traces of gamma matrices 



and other manipulations were done using program TRACER [|16| written 



for MATHEMATICA. Dimensional regularization has been used and the 75 



matrix has been treated according to the 't Hooft-Veltman scheme |17 
The result for the twist-3 contribution has the general form 

Tr[^7 5 lWp)] = 

A^f (-p 2 f\og{-p 2 /^ 2 ) + B^< qq > 2 log(V//i 2 ) 



-C- 4 < fg 3 GGG > log(-p 2 /V) 



n , n, / 2 / 2^ a S mo < qq >' 
D + Flog(-p j\i ^ 1 



7T — p 2 
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A 


B 


C 


D 


F 


G 


H 


d 


l/(3 4 2 7 ) 


8/27 


l/(3 2 2 6 ) 


85/162 


-95/648 


10/81 


2/9 


u 


l/(3 4 2 6 ) 


4/9 


l/(3 2 2 5 ) 


-847/2592 


-13/324 


2/81 


-1/27 



Table 1: Numerical coefficients corresponding to the sum rule eq. [16]. The 
upper line gives the values for the twist-3 operator involving d quarks, the 
lower line for the corresponding operator with u-quarks 



< qq> 2 < g 2 G 2 > < qq > 2 m\ 
+(jr : ttzt, V ti - 



-p 2 ) 2 (—p 



,2\2 



S[aP{n]Pu} (16) 



where the numerical coefficients are given in Table [I]. Note that the coefficient 
in front of the gluon condensate (GG) vanishes. 

There are two types of unwanted terms at the phenomenological sides 
of the sum rules. The first type are the terms containing a {m 2 N — p 2 )^ 1 
singularity. The second are pure continuum contributions which correspond 
to contributions from various higher resonances etc. To deal with these 
terms we first multiply both sides of the sum rules by (rn 2 N — p 2 ) and then 
perform the Borel transform and a suitable continuum subtraction. The 
multiplication has a twofold effect: (a) it lowers the singularity of the first 
term and (b) removes the m 2 N — p 2 singularities from all other terms. In that 
way the only low lying singularity is that of the first term. The singularities 
at higher values of — p 2 are suppressed by the Borel transformation. The 
transformed sum rule for twist-3 contribution looks as follows: 

-lQd^Xj\ G m%e-</ M2 

= 3!M 12 i(^|£ 4 - 4£ 5 ) + M\B + C)(E 2 - ^E 1 ) + M 2 m 2 N D + 
M4 F(^ MM 2 //, 2 ) - Ey) +M 2 (G + H)(^ + 1) (17) 



where 

n— 1 i , 

E n = l- e (~ so/M2 ) Tj (so/M 2 ) (18) 



k=0 

Here we absorbed the condensates, as and tt factors in the definition of the 
coefficients A,B, . . . which then are given by A = as A/it 5 , B = as < qq > 2 
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B/tt, C =< fg 3 GGG > C/ir 4 , D = a s m 2 < qq > 2 D/tt, F = a s m 2 < qq > 2 
F/ir, G =< g 2 G 2 >< qq > 2 G, H = m 4 < qq > 2 H. 

One point requires explanation here. The term associated with the dimension- 
8 condensate in the sum rule for the matrix element of the twist-3 operator 
contains arbitrary scale /i even after Borel transformation. This dependence 
may be traced back to the occurence of terms proportional to log(— p 2 / fj, 2 ) /p 2 
in the original sum rule. They arise because of the mixing of the three-point 
correlation function (|5|) with two-point correlation functions of the nucleon 
currents rji(x) and fja(0) with operators produced by the contraction of Or{y) 
with fjaiO) and ni(x), respectively |Tj|. The natural choice of the scale is // ~ 



1 GeV which corresponds to the typical values of momenta in the intermedi- 
ate region where the sum rule equality is expected to hold. Fortunately the 
numerical coefficient associated with this term is small. 

The condensates have been normalized at 1 GeV. Their values: < qq >= 
(-0.257GeV) 3 , < a s /nGG >= 0.012GeV 4 , m 2 =< qgaGq > / < qq >= 
0.65GeV 2 and < fg 3 GGG >= 0.046GeV 6 correspond to the standard ITEP 
values rescaled to the normalization point fi 2 ~ m 2 N ~ 1 GeV 2 . The strong 
coupling at 1 GeV is taken to be as = 0.37 (A = 150) MeV. The continuum 
threshold is taken at so = (1.5GeV) 2 roughly corresponding to the Roper 
resonance position. 

The coupling constants Xj and Xg may be calculated from the two point 



sum rule derived in 13 



2{2nYm 2 N X 1 X G e' m -/ M2 = -^M 8 E 4 

5 7T 

+ 1 < g 2 GG > M 4 E 2 - ^(2tt) 4 < qq > 2 M 2 E X + l(2n) 4 m 2 < qq > 2 
2 3 7T 3 

(19) 

However, it is more convenient (and usually more accurate) to divide simply 
the sum rule ( |T7p by the sum rule (|T^). The coupling constants at the 
phenomenological side cancel out. The quotient sum rule for the matrix 
element of twist-3 operator has been plotted in Fig. [I]. This figure shows 
the singlet (S) and nonsinglet (NS) part of S 2 \ S 2 \S) = dP\u) + S 2 \d), 
dp) [NS) = (u) — (d) . The corresponding matrix elements d^ (proton) 
and d^ (neutron) are shown in Fig. 0. For comparison the corresponding 
sum rules obtained from the the analysis in || which employed Ioffe currents 
only are shown. The corresponding coupling constant A 2 is determined from 



9 



0.2 



0.15 



s 1/2 = 1.5 GeV 
/j 2 = 1 GeV 2 



0.1 

& 0.05 



CM 




Figure 1: Stability plot of the sum rule eq. (0). The full line corresponds 
to d^(S), the dotted line to d^ 2 \NS). For comparison the results of the 
analysis in M are also shown. The dashed line corresponds to d BBK (S) the 
space-dotted line to d BBK (N S). 
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c 
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1.6 



Figure 2: Dependence of S 2 ' (proton) and S 2 ' (neutron) on the Borel pa- 
rameter M 2 . The full line correponds to the S 2 ^ 1 the dotted line to S 2 \ The 

(2) 

dashed and space-dotted lines correspond to the results for d^ K (proton) and 
dftft K (neutron) obtained from the sum rules given in ||. 
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the additional two point sum rule 



2(2n) 4 Xje-<^ M2 = M 6 £ 3 

+ i < g 2 GG > M 2 E 1 + ^(2tt) 4 < qq > 2 (20) 

which is just the standard sum rule considered by Ioffe. Instead of using a 
fixed value for X 2 we divided the sum rule obtained in || by the sum rule 
( |20|) which improved the stability considerably. The final values of the matrix 
elements may then be estimated from the figures by taking M 2 ph m 2 N : 

S 2 \S) = -0.068 ± 0.03 d (2) (NS) = 0.078 ± 0.03 
d {2) (proton) = -0.006 ± 0.003 

d {2) (neutron) = -0.03 ± 0.01 (21) 

Errors are due to the M 2 dependence of the sum rule, see Figure |l] and || 
resp. The uncertainty introduced by the hypothesis of factorization of higher- 
dimension condensates is discussed below. These values are to be compared 
with those obtained from the sum rules given in M: 

d ( zl K (S) = -0.05 ± 0.03 d<£l K (NS) = 0-065 ± 0.03 
dBBKiproton) = -0.003 ± 0.003 
d {2 l K (neutron) = -0.025 ± 0.01 



In a recent paper |19| Ross and Roberts have estimated contributions to 
the sum rule ( |T5| ) arising from chiral corrections to the phenomenological 
spectral density. However it is not clear how to take into account the cor- 
responding chiral contributions on the theoretical side of (|ToD, see e.g. PD| ]. 
When the necessary "theoretical error" is taken into account, their result 
agrees with that of ||. As this error is unavoidably large also in the present 
case we have not applied the Ross and Roberts analysis to our results. 

As in H and [18| our sum rule is dominated by the contribution from 



the highest dimension operators considered i.e., those with coefficients de- 



termined to the leading accuracy by tree diagrams. As discussed in [18] this 



could in principal signal a breakdown of the OPE for the correlator in ques- 



tion which may not allow to make a reliable numerical estimate. In |18| this 



potentially important problem was settled by an estimate of the contribution 
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from operators of yet higher dimension 10. It was found that the contribution 
of operators of dimension 10 is likely much smaller than those of dimension 
8, which dominated the sum rule in [[TJJ. This is a situation very typical 
for QCD sum rules. Low-dimension operators are usually numerically unim- 
portant because each extra loop in the definition of the coefficient in front 
of an operator of lower dimension brings in a small factor ~ (27T 2 ) -1 . On 
the other hand the higher the dimension the stronger is the suppression due 
to the Borel transfromation. The dominant graphs are therefore usually the 
lowest dimension tree graphs. In the present case these are the dimension 10 
contributions. We expect that as usually still higher dimension contributions 
starting with dimenision 12 will be small due to the Borel transformation. 
Unfortunately a quantitative test of these arguments is difficult as very little 
is known about VEV of operators of dimension 12. 

Note also that our calculation involves factorization of the condensate 
of high dimension. In this respect it is also very similar to the calculation 
presented in ref . || . The factorization of the higher dimensional condensates 
is a standard technique used in the QCD sum rules. It is necessary to reduce 
the number of unknown parameters. The physical assumption hidden behind 
this procedure is the absence of certain higher order correlations in the QCD 
vacuum. The validity of this assumption has been studied for several simpler 
condensates pT| . Yet, the knowledge of condensates of dimension 10 is scarce 
because such condensates rarely appear in calculations. One can even turn 
this argument around and argue that the good agreement with the previous 
calculation [Q supports the validity of factorization hypothesis. The BBK 
result may be treated as a kind of cross check because it relies on a different 
condensate for the dominant contribution, namely < qq > 2 ??1q instead of 
< qq > 2 < G 2 > and < qq > 2 mj. The other test is the two point sum rule 



for the nucleon mass with two T)g currents [13|]. This calculation involves 
condensates of the dimension 10 and leads to reasonable results. 

The agreement between the two calculations is a highly non-trivial check 
of the QCD sum rules method. The main difference is that in || gluons were 
generated entirely by perturbation theory while in our case they are already 
present in the interpolating current. As a consequence our correlator (|5[) has 
milder singularities than that analysed in ||. In principle this should be an 
advantage of our approach but as the current r] G (x) has two units higher 
dimension than that of the f]i{x), the accuracy is diminished by the need to 
consider higher dimensional condensates. Given an approximate character 



13 



of the method it is by no means obvious that both approaches give the same 
results. It seems to us that our results not only give a strong argument in 
favour of the numerical estimates presented in fl2"T| ) and (|2"2]), but also give a 



strong support to the credibility of QCD sum rules method. However, neither 
type of sum rule can be expected to yield very accurate results, because of 
the slow convergence of the operator product expansion. 

We think that taking all the uncertainties into account the following three 
qualitative statements can be made: 

• \d^\ and \d^\ are smaller than 0.03. The corrections from the twist-3 
operator to the Bjorken sum rule and Ellis- Jaffe sum rule are therefore 
of order 0.013m^/<5 2 which e.g. for / g\(x)dx amounts to a correction 
of ~ 5% at Q 2 ~ 2GeV 2 . Note that for a definite answer also and 

have to be taken into account. 

• There is a strong isospin asymmetry between d)?> and d^\ \d^\ > 
\dffl\. This isospin asymmetry should manifest itself in a number of 
spin phenomena [T(| . 



• the sign of df^ is negative (in contradiction to the bag model prediction 
0). d^ is negative. 
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